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THE EVALUATION OF QUANTUM INTEGRALS 

By Edwin C. Kemble 

Jefferson Physical Laboratory, Harvard University 

Communicated by E. H. Hall, April 22, 1921 

The application of the Wilson-Sommerfeld quantum conditions to a 
conditionally periodic system with orthogonal coordinates involves the 
evaluation of an integral of the type 

J = f <JWdq. (1) 

The integral is to be extended over a. complete cycle of values of q, 
which oscillates between two roots of f(q). The sign of the radical is to 
be the same as that of dq, so that if a and b denote the roots of f(q), the 
integral can be written 

J = 2 ^Wdq (2) 

J a 

Iif(q) is a polynomial of the second degree in either q or 1/q the integral 
can be cleanly evaluated. Otherwise, approximations are generally 
necessary. If J(q) can be expressed in the form 

/(«) = <p(q) + « t(q) 
where <p(q) is quadratic in q or 1/q, a is constant, and a\p(q) is small, a 
natural method of procedure is to try to develop J into a power series in 
a. Thus 

/(0) and 7'(0) are easily evaluated, but unfortunately the higher deriva- 
tives of J with respect to a cannot be calculated by the usual methods be- 
cause the higher derivatives of V(/g) with respect to a become infinite 
at q = a and q = b. Hence this method is useful only when the higher 
order terms are negligible. 1 

Another method of attack employed by F. Tank 2 and accepted as 
valid by others 3 turn- it on close examination to be faulty. Tank de : 
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velops f(q) into a power series about its maximum point, q = d. Let 
£ = q — d, and let # denote the maximum value of /(g). Then f(q) can 
be thrown into the form 

f(q) =H-a?- {$?+ v? + 8£* + . . . ). (4) 

Introducing the symbols 

V(£) = «?, 

AV = /5f» + 7f* + «{' + ..... 
we obtain 



V/(a) = Vfl-F^l- 



AF ,— — . 1 AF ' 1 (AF) 2 



H-V 2 V-ff-V 8 V(i7-F> 

Tank in effect integrates this series term by term between the limits 

H *A 

— (H/a) and +(H/a) and identifies twice the sum of the series so obtained 

with J. This procedure is wrong, since the correct limits of integration 
for J are £ = a — d and £ = b — d. Moreover, the expansion is not usually 
convergent throughout the interval of integration, so that it is not possible 
to correct Tank's work by altering the limits. 

Another scheme of series development may be suggested, which avoids 
the above difficulty. Let the quantity u be defined by the equation 

H = V h - /(g) = £ ^^+Pi+yp + sf+.... (5) 

The sign of u is to be the same as that of £. The integral J can now be 
thrown into the form 



J = 2UH-u^du (6) 

I du 

u = — VF 
Let us assume that d%/du can be developed into a power series in u. Thus 

00 



du = 2 ay (7) 



r=*0 



Let f+VU 

£T T = 2\ » T < H ~ u *du. (8) 

J-Vh- 

To evaluate # T we introduce the variable of integration 6 defined by the 
relation u = tf !i sin 8. 
Then (8) becomes 

K T = H 2 \ sin T 6 cos*0 d». (9) 
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It is easy to show from (9) by the application of well-known formulae that 
K T vanishes for odd values of r; that K is tH; and that for even values 
of r greater than zero K r is given by the equation 

K T = 2ir H 2 <Kt) (10) 

where 

1.3.5.... (t-1) (11) 

« {T) = 2.4.6.... (r+a)" 
Now (6) becomes m 

J = tH [ a + 2 ^ a 2r^ (2r)H T ] . (12) 

We proceed to the evaluation of the coefficients a , oi, o 2 , .... As d£/d« 
may be regarded as a function of either £ or u, let 

df/d* = f(u) = x ft). 

Then (T) 

a r = i ^(0). (13) 

Since » and £ vanish together, the derivatives of ^ at the point u = 
can be calculated from the derivatives of X at the point £ = 0. The 
following method of procedure is perhaps the simplest. Let 

*<*> = ^p- = « + « + 7f 2 +... (14) 

Then « = £ - A and - =— - Ep^-i . 



Let w-2«(«)+f«'({)=---^i±*.. 



(14') 



Then x©=?= ^. (15) 



Differentiating, 

* '(«) = X' (I) X ft) 



*■(-)(«) -xft) |{^VJ 



(16) 



In some cases these successive derivatives of ^ are simple and easily 
calculated functions of £. In others the successive derivatives rapidly 
become complicated. If f(q) is given as a power series in £ the process 
of differentiation can be performed conveniently as follows : 
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f («) = — , 



**(«) = 



f'"(u)=~\W + 2 W '^^L 



4Vt> / » _ 4»'w' 2W6* 6»(w') ! 
t» 3 I w ww 2 

Uvv'W 4 (»')*»' 14 wV 



50wV) ! 40*t ! wV' _ 4t>W ' __ 60 »«(w')' 



(17) 



At the point 

«=J=0, 

at = 2»=2«, w' = 3»' = 3/S 



M) >'=4j ) ''=8y 
w (^ = (n + 2)» w = 

Hence 



(18) 



*(0) = 



v; 



^ --£{;*-. 



(V 2 I 4 
" «'A 1 16 a* 2 a 2 

i with 



+ 



4.1 



(19) 



a'/s I 16 a 8 2 a 2 ' a ' a 7 / 

Combining the above equations with (13), (11), and (12), we obtain 
J = r i> (o)H + 2 x ? (2) — -F 2 + 24 ,, (4) '^p-H 3 +... 

+ — +~--e Lff 2 + I 

a a 7 J I 

This result is similar in form to Tank's and the coefficients of the two 
lowest powers of H are the same as in his series, but the coefficient of H* 
is different. 

In order that the result may be valid, it is necessary that the series (4) 
converge and represent f(q) throughout the interval of integration, and 
that the derivatives of ^ shall all be finite throughout a circle whose radius 
is greater than -\~H drawn about the origin on the complex u plane. From 
(17) it is evident that the derivatives of \p will be finite up to the point 
where w vanishes. It follows from (14') that/' (#)/(<? — d) must have no 
zeros on that part of the complex q plane which is mapped on the above- 
mentioned circle on the u plane. It is necessary, in particular, that 
f(q) shall vanish not more than once for real values of q between a and b. 
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It does not seem worth while to attempt an exact discussion of the 
boundaries of the region on the complex q plane from which the zeros of 
/'(?)/(9~^) must be excluded, but we can say qualitatively that there is 
little chance that the series (7) will not converge throughout the interval 
of integration if the series 

w = 2« + 30£ + 4y£ 2 + 5S£ 3 + . . . (21) 

converges rapidly. 

In the applications of this analysis an expression for H as a function of 
J will generally be desired. The power series may be reversed to ad- 
vantage by the following scheme which resembles that of Lagrange. Let 
H = F(J). Then by Taylor's theorem 

H = JF'(0) + ~F"(0) + J -F (0) + .... {22) 

Let b n denote the coefficient of H n in (20). F'(0), F"(0), etc., can be 
calculated in terms of the b's. Let us first compute the derivatives of 
F in terms of H. Let 

co 

P (H) = 2^« H "~ 1= W 

n — \ 

Then 

dH 1 

Tj = F 'U) = W) - (23) 

Differentiating again, 



F'U) 



dH 



f-\ EE _=£.' 

\PJ'dJ ~p* 
V<>,< r) = ± ( _A d E _ P" i 3(»') 



(24) 



Since / vanishes when H does, the values of F'(0), F"(0), etc., are ob- 
tained by setting H equal to zero in the right-hand members of (23) and 
(24). Thus 

F'(0) = l/h; F"(0) = -2h/bl; 

F'"(0) = - 6b 3 /bj + Ylb\lb\; 



(25) 



If these formulas for F'(0), F"(0), etc., are evaluated in terms of the ex- 
pressions for the b's given by (20), equation (22) becomes 

r 8x 2 a[_' 4«J ' 32*- 3 a A J a 
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2 a 4 a 2 32 a 3 J + - • • (26) 

As a check on the series development here suggested, the writer has 
derived Sommerf eld's formula, 



W' + ^«.-.^c[;fe-.] 



(27) 



by means of the series (12). The derivatives of ^ when evaluated by- 
equations (16) are in this case particularly simple and as a result the series 
can be summed. 

It is of interest to note that the method of development in series here 
suggested is applicable to a variety of problems. It may be used to 
evaluate indefinite as well as definite integrals. 

One simple application is in the determination of the periods of os- 
cillation of a vibrating system. Iu the case of a conditionally periodic 
system with orthogonal coordinates, the periods are given by expressions 
of the form 

dq 

where a and b are again roots of f(q). Introducing the quantities H, |, u 
defined as in the preceding work, we obtain 

S di i du 
, du. (29) 

_^H-u* 

If d£/du is developed into a power series of the type (7), we obtain 

00 
T = 2 ^r 

where 



^ *? = 2 fe (SO 




-Vi? 

The coefficients K T vanish when t is odd, and the final expression for the 
period is 

00 

T = fcr [ a + 22^r (r+1) V (2t) ff J . (30) 

T = l 

1 This is in effect the method used by Sommerfeld, Atombau und Speklrattinien, 
2nd edition, Braunschweig, 1921, pp. 476-182. 

a F. Tank, Mitteilungen Physik. Ges. Zurich, 1919, p. 87. 
« Cf. G. Hettner, Zs. Physik., 1, 4, 1920 (350). 



